Flashing annihilation term of a logistic kinetic as a mechanism leading to Pareto 

distributions. 
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It is shown analytically that the flashing annihilation term of a Verhulst kinetic leads to the 
power-law distribution in the stationary state. For the frequency of switching slower than twice 
the free growth rate this provides the quasideterministic source of a Levy noises at the macroscopic 
level. 



I. INTRODUCTION 

1// noise [1-6] and power-law distributions [6-9] are 
widely observed signatures of anomalous behavior [7,10]. 
Such phenomena are related to the scale invariance of 
considered systems and appear in different domains of 
natural and social sciences. Their universal character is 
partially explained on the statistical manner as a conse- 
quence of the central limit theorem for scaleless processes, 
where the Gaussian distribution is replaced by appropri- 
ate a-stable distribution [7,10]. This corresponds to the 
use of the Levy noises to describe fluctuations or frac- 
tional dynamics at the level of Fokker-Planck equation 
[11,12]. Anomalous behavior can also be obtained for 
nonlinear systems driven by (conventional) multiplica- 
tive Gaussian white noise [8,6], particularly for Lotka- 
Volterra (or Malthus- Verhulst) model [8]. In any case 
processes exhibiting fat-tailed distributions are subjected 
to large fluctuation, so within certain stochastic descrip- 
tion it is difficult to distinguish to which extent such 
behavior is an intrinsic property of the kinetics and to 
which extent it is simply the reflection of the presence 
(and specific properties) of noise. The source of the pos- 
sibly intense noise is known for thermodynamical critical 
systems. The search for explanation of scale invariant 
phenomena in generic systems is still unfinished [2-6,8]. 
Particularly, it is suggested in Ref. [2], see also Ref. [5], 
that the 1// noise in a membrane channel is produced 
by random switching between both conducting states, 
rather than by inherent properties of ion transport. Simi- 
larly, the importance of activation-deactivation processes 
is addressed in Ref. [3] . 

In the present paper we are going to show analyt- 
ically that the asynchronic switching between generic 
pure Malthusian and Verhulst-type kinetics leads to the 
Pareto distribution in the stationary state. The result 
seems interesting for few reasons. The logistic equation 
is used, which is the basic model of evolution for social 
sciences. The power-law distribution was identified by 
Pareto to describe the social statistics, namely the dis- 
tribution of (large) wealths. The model is quasideter- 
ministic, so the asymptotic power-law behavior has the 
kinetic origin only, coming as a result of balance between 
Malthusian growth and Verhulstian saturation. Finally, 
our result supports the general opinion of Ref. [2] about 
possible sources of anomalous behavior. 



II. MODEL 

Power-law behavior may be considered as an interme- 
diate one between divergent and bounded kinetics. Let 
us note that the Malthus equation, 



x t 



ax, 



(1) 



where a is a positive difference between birth and death 
rate, leads to the exponential growth of population. Tak- 
ing into account the Verhulst competition coefficient bx^ , 
\i > (and equal unity for the "true" Verhulst model), 
depending on the actual size of population, 
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bx» +1 



one obtains, instead of unbounded trajectories, 
monotonic relaxation to the stationary value 

a; s t = {a/by, v = 1/fi, 

with the relaxation time 

T = 1/fia. 



(2) 
the 
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(4) 



So, if the annihilation term of Eq. (2) is absent (b = 0) 
the kinetic is divergent and as long as b > the kinetic is 
bounded. Thus one expects that the "flashing" annihi- 
lation term, temporarily switched on and off, may result 
in a power-law distribution for the asymptotic state. Let 
us simply assume a (Markovian) binary and asynchronic 
character of a switching process 



b(t)=b[l + al t ], 



(5) 



where I t — (— 1)"^*, where Nt is a Poisson process with 
parameter A and a = +1 or — 1 with probability 1/2. 
It means that the evolution consist of periods of active 
annihilation (with a coefficient 2b) separated by periods 
of a free growth. The length of the periods is random 
with the average equal A -1 . 



III. CRITICAL cr 2 — DEPENDENCE OF 
MOMENTS 

The Bernoulli equation (2) driven by dichotomous 
Markov process has been already considered [13,14] . Par- 
ticularly, in Ref. [14] the transient behavior of (2) with 
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nonlinear coupling (5), \a\ < 1, has been examined and in 
the limit t —>■ oo the following formula for the stationary 
moments 

(^^^(^f^l + c^ 2 ) (6) 



where 



(5 = X/fm- 1/2, 



(7) 



was obtained. 

The value a 2 = 1 is critical with respect to the con- 
vergence of the hypergeometric series (6) (branch point 
of the function) and for the critical exponent 



1 + S - sv - 1/2 = (A/a - s)v < 0, 



(8) 



the r.h.s. of Eq. (6) diverges at a 2 = 1. The condition 
(8) clearly shows the existence of fat-tail 



P(x) ~ ar 1 "*/" 



(9) 



for large x. 



IV. STATIONARY PROBABILITY DENSITY 

The stationary probability density is formally given by 
the inverse Mcllin transformation (with respect to s) of 
Eq. (6). The more convenient way of computation is 
to apply the so-called quadratic transformation of the 
hypergeometric function [15] 



{x% t = x s st (l + a)-» s 2 



Frfa^+S; 1+25;^-) (10) 



and then use of Euler's integral representation 
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5(7,7) 
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l + a 



(ii) 

where 7 = 1/2 + 5 and B is the Beta function. Introduc- 
ing a new variable £ = [6(1 + cr — 2ert)/a]~" wc obtain 
finally 



where 

P(x) = N-'x-^ 1 [a 2 - (1 - a^/6) 2 ] ^ 1/2 , 
and x G (x/, x r ), 



(12) 



(13) 
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i6(l + M). 







(14) 

The normalization constant in Eq. (13) is equal ./V = 
S(7, 7)(2|cr|) 2l5 6/(^a). Note that in the case of interest, 
I cr I = 1, the right boundary x r — > 00 and 

P(x) oc x- 1 - x ' a {2 - ax-^/bf- 1 ' 2 , (15) 

where x € ((a/26)", 00), in agreement with Eq. (9). 



V. NOISE-INDUCED TRANSITIONS 

Let us remind that in the present work we consider 
Eqs. (2), (5) with \a\ — 1 as a noiseless evolution consist- 
ing of two essentially different (and randomly switched) 
modes of deterministic kinetics. On the other hand the 
same model, especially for \a\ <C 1, may be treated as 
a stochastic kinetic with nonlinearly coupled paramet- 
ric noise, when the fluctuation of a coefficient 6 are de- 
scribed by a zero mean dichotomous color noise £ t = 60-/4 , 
(£t£o) = 6 2 er 2 e~ 2At , of a correlation-time r and intensity 
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i2 2 

ban 



(16) 



From this point of view it is worth to analyze the r and D 
dependence within the context of so-called noise-induced 
phase transitions [13]. The case a 2 = 1 correspond to 
the critical line b~ 2 D/r = 1 in Fig. 1 and to the "true" 
transition, which is reflected by singular cr 2 -dependence 
of observables, Eq. (6), and related power-laws. Be- 
low this line the kinetic is bounded and consequently 
the support of stationary distribution is finite, Eq. (14). 
Above the line the system is unstable and x t rapidly ap- 
proaches infinity (after finite time). At the critical line 
the stationary states are the semiaxis distributions (lower 
curves in Fig. 1), Eq. (15), exhibiting fat-tails with in- 
dices a — A/a. The dependence on the correlation-time 
is also interesting. If the frequency of switching A is small 
(long correlations), 5 — 1/2 = A//za — 1 < 0, both bound- 
aries xi and x r — which are the stationary values (3) 
of a deterministic kinetics (2) with coefficients 6(1 + |ct|) 
and 6(1 — |cr|) respectively — are attractive, the proba- 
bility density (13) has a minimum and it is convex down. 
In contrast for A greater than the relaxation rate /xa, 
Eq. (4), the boundaries become repulsive. P(x) is then 
unimodal. The A = fia [or r = l/(2/za)] is a critical value 
with respect to the noise-induced transition in probabil- 
ity density. 



VI. EXACT PARETO DISTRIBUTION 

The case A = fia, see the vertical line in Fig. 1, is ex- 
ceptional for other reasons. The stationary distribution 
is then purely power, P(x) oc and neither diverges 

nor vanishes at the boundaries (14). And for \a\ = 1, at 
the critical point common for both lines, it becomes the 
exact Pareto distribution 

P(x) = (afi/2b)x- 1 ->*, x > (a/26)". (17) 

Moreover the case A = fia was the nontrivial one for 
which the time-dependent transient moments were found 
in the closed analytical form [14] 

(xl) = e-» at {x s (t;0)+x s (t;2b)}/2 + ^- b 

x ij^[x s -^t;2b)-x s -^t;0)}, if » ± s 
|slog[a;(t;0)/a;(i;26)], if \i = s { ' 
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where 

x(t; 6) = [xq "e-"°* + 6(1 - e-"°*)/a] (19) 

is the solution of Eq. (2) for a constant a and 6 param- 
eters. Using the asymptotic forms x(t; 0) = xoe at and 
x(t;2b) — > (a/26)", one concludes that the higher order, 
s > yu, transient moments diverge exponentially 

(x t s ) oc e («-M)«t (20) 
and for the marginal case s = n (= A/a) 

2 

(a;?)- const (21) 

grows linearly for long times. Note that according to 
(15) or (9) the value of the ratio a = A/a (frequency 
of switching by the rate of a free growth) specifies the 
order of the lowest divergent moment. In the "double- 
critical" case of exact Pareto distribution (17) it is simul- 
taneously equal to the value of a "degree of nonlinearity" 
/U. Thus, for the most important true Verhulst or so 
called Stratonovich (with /i = 2, corresponding to the 
quartic potential) model it means that it is the mean 
value or the variance, respectively, that grows linearly 
with time, Eq. (21). 

VII. CONCLUSIONS 

The anomalous systems exhibit large fluctuations 
which arc difficult to explain dynamically within stan- 
dard theories, for instance, as a result of an additive 
Gaussian diffusion. Thus such a behavior is frequently 
joined with the presence of a "large" driving noise, e.g., 
a Levy noise or a multiplicative Gaussian noise. On the 
other hand it was suggested that certain aspects of such 
behavior are rather related to some hierarchical struc- 
ture of kinetic processes at the deterministic level [2,3,5]. 
This idea was exploited in the present work for exactly 



solvable generic model (2), (5). The general behavior of 
the system depends on two things, on the mutual relation 
between A and a (or /j,a) and on the value of a 2 , namely 
whether a 1 < 1 or a 2 = 1. 

Let a 2 < 1 (i.e., below the critical line in Fig. 1). If 
frequency A is smaller than the rate of relaxation pro- 
cesses [ia the x t "hardly" follows b(t), Eq. (5), switching 
between some values located close to xi or x r , Eq. (14), 
respectively (right part of plot). Successive ("rare" and 
"short" ) moves up and down are associated with values 
6(1 — |ct|) and 6(1 + |cr|) in Eq. (2), respectively, and in 
principle look like different deterministic processes. In 
the opposite case (left part of plot), A > /j,a, the evolu- 
tion looks more li ke th e Verhulst kinetics with the av- 
eraged value 6 = b(t) of annihilation coefficient. The 
distribution of x is unimodal, however still with a rather 
wide maximum, located somewhere in a middle of the 
support. The value A = \ia is critical for so-called noise 
induced transition. This transition is formally controlled 
by correlation-time of the noise. 

The case of a 2 = 1 is of particular interest for the 
present work and may be formally considered as the 
"true" transition, which is indicated by singular depen- 
dence of moments and related power-laws. This transi- 
tion is controlled by the noise intensity. The asymptotic 
distribution of x is given by Eq. (15). P(x) has a fat-tail 
with the index a = A/a. The distributions with a < 2 be- 
long — via central limit theorem — to the basin of attrac- 
tion of appropriate Levy ev-stable distribution. Thus the 
kinetic (2), (5) with randomly flashing annihilation term 
may be treated as a source of a Levy noise at the macro- 
scopic level, if the frequency of switching is small enough, 
A < 2a, or, cquivalcntly, if the correlation-time is suffi- 
ciently long t > (4a) _1 . For the Stratonovich model, 
which is presented in Fig. (1), this corresponds precisely 
to positions on the line right from the critical point [with 
t = (2/ia) -1 ]. If the frequency is large, A > 2a, or the 
correlations are short, P(x) still has a fat-tail, however 
it corresponds to the Gaussian distribution via central 
limit theorem. 
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